Asymptotic Safety of Gravity Coupled to Matter 
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Nonperturbative treatments of the UV limit of pure gravity suggest that it admits a stable fixed 
point with positive Newton's constant and cosmological constant. We prove that this result is 
stable under the addition of a scalar field with a generic potential and nonminimal couplings to 
the scalar curvature. There is a fixed point where the mass and all nonminimal scalar interactions 
vanish while the gravitational couplings have values which are almost identical to the pure gravity 
case. We discuss the linearized flow around this fixed point and find that the critical surface is 
four-dimensional. In the presence of other, arbitrary, massless minimally coupled matter fields, the 
existence of the fixed point, the sign of the cosmological constant and the dimension of the critical 
surface depend on the type and number of fields. In particular, for some matter content, there 
exist polynomial asymptotically free scalar potentials, thus providing a solution to the well-known 
problem of triviality. 
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I. INTRODUCTION 

The failure of perturbative approaches to quantum 
gravity does not necessarily imply that quantum gravity 
does not exist as a field theory. There is still in principle 
the possibility that the theory could be "nonperturba- 
tively quantized". To understand what this means, one 
has to look at the Renormalization Group (RG), i. e. the 
flow of the coupling constants gi{k) as a certain external 
momentum parameter k is changed [22]. It is customary 
to take k as unit of mass; if gi has dimension di in units 
of mass, we define dimensionlcss couplings gi — gik~'^^ . 
The RG flow is then given by the integral curves of a 
vectorfield (3 in the space of all couplings, whose compo- 
nents Pi{g) ~ dtgi (with t = Infc) are the beta functions. 
A Fixed Point (FP) is a point in the space of all cou- 
plings where 



A(<?*) = 



(1) 



For example, in ordinary quantum field theories in 
Minkowski space, the point where all couplings vanish 
is always a FP (called the Gaussian FP), because a free 
field theory does not have quantum corrections. 

Suppose that the theory admits a FP. One defines the 
critical surface to be the locus of points that, under the 
RG evolution, are attracted towards the FP for t ^ oo. 
Starting from any point on the critical surface, the UV 
limit can be taken safely, because the couplings, and as a 
consequence the physical reaction rates, will be drawn to- 
wards the FP and hence remain finite [23] On the other 
hand, if one starts from a point not belonging to the 
critical surface, the RG evolution will generally lead to 
divergences. If the critical surface has finite dimension c. 
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the theory will be predictive, because only c — 1 param- 
eters will be left undetermined and will have to be fixed 
by experiment at a given energy scale (the last remain- 
ing parameter being the scale itself). The special case 
when the FP with a finite dimensional critical surface is 
the Gaussian FP is equivalent to the usual perturbative 
notion of rcnormalizability and asymptotic freedom. 

This scenario for nonperturbative rcnormalizability 
has been discussed specifically in a gravitational context 
in [1] where this good property was called "asymptotic 
safety". At the time some encouraging results were ob- 
tained by studying gravity in 2 -f e dimensions [1, 2], but 
the program soon came to a halt essentially for want of 
technical tools. It now appears that the right tool to 
tackle this problem is the Exact RG Equation (ERGE), 
in one of several guises that have appeared in the litera- 
ture in the last decade [3-5]. 

The ERGE is a differential equation that determines 
the RG fiow of the action. It can be viewed as a set 
of infinitely many first order differential equations for 
infinitely many variables (the coupling constants), and 
therefore cannot generally be solved in practice. A 
method that is commonly used to calculate nonpertur- 
bative beta functions is to make a physically motivated 
Ansatz for the running effective action, typically contain- 
ing a finite number of parameters, and insert this Ansatz 
into the ERGE. 

The ERGE, in the specific form discussed in [5], has 
been applied to Einstein's theory in [6, 7], where the beta 
functions for Newton's constant and for the cosmological 
constant were derived. It was later realized that these 
beta functions actually admit a nontrivial UV-attractive 
FP [8]. The properties of this FP were further discussed 
in greater detail in [9]. A particularly important issue 
is to prove that the FP is not an artifact of the Ansatz 
but is a genuine property of gravitation. Several facts 
seem to indicate that the FP is quite robust. It has been 
shown to exists in four spacetime dimensions for many 
different shapes of the cutoff, whereas in other dimen- 
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sions it only exists for certain cutoffs but not others. Its 
properties have been shown to be only weakly dependent 
on the shape of the cutoff, indicating that the trunca- 
tion is self-consistent [10]. It has been found in certain 
dimensionally- reduced versions of the theory [11]. It is 
also remarkable that Newton's constant always turns out 
to be positive, a fact that could not in any way be guar- 
anteed by the general form of the equations. 

The most important test, however, is the stability of 
the FP against the addition of new couplings. Each time 
we consider a new coupling, whether remaining in the 
context of pure gravity or if we introduce matter fields, 
a new beta function has to vanish and therefore a new 
constraint has to be satisfied by the set of all couplings at 
the FP. It is therefore nontrivial that the FP still exists 
when we take into account additional couplings. 

In the context of pure gravity, an important progress 
was made in [12], where it was shown that the addition 
of a term quadratic in curvature does not spoil the ex- 
istence of the FP. In fact it turned out that the values 
of the cosmological and Newton's constants at the FP 
are almost unaffected by the new interaction, while the 
new coupling constant is quite small at the FP. This is 
far from being conclusive evidence, but it is nevertheless 
an important result, especially in view of the fact that 
another FP that was present in the truncation with only 
two couplings — the Gaussian FP — does not exist in the 
three coupling truncation. 

As far as matter is concerned, wc have recently consid- 
ered the effect of minimally coupled, masslcss quantum 
fields of arbitrary spin [13]. The only couplings taken into 
account were the cosmological and Newton's constant, 
since the coefficients of the matter kinetic terms can be 
normalized to their standard values by field rescalings. It 
was shown that the existence of the FP, the values of the 
cosmological constant and Newton's constant at the FP 
and the dimension of the critical surface all depend on the 
type and number of fields present. Altogether, the exis- 
tence and attractiveness of the FP puts some constraints 
on the number of matter fields that are present. 

In this paper we continue the analysis of coupled grav- 
ity and matter systems and we begin to address the issue 
of matter couplings. There are many different couplings 
that are necessary to construct realistic theories of the 
world, and we cannot possibly take them all into account, 
so as a first step we shall consider the simplest example, 
that of a self-interacting scalar field. Aside from its role 
as a model for the Higgs field in unified theories, a scalar 
field (the dilaton) appears in many popular theories of 
gravity. It can therefore sometimes be regarded as part 
of the gravitational sector, rather than the matter sec- 
tor. This makes its properties especially interesting in a 
gravitational context. In this paper we shall not make 
any assumption about the physical interpretation of the 
scalar field. 



The class of actions that we consider is 

T[g, <t>\^ j d^x-Vs {v{cf') - F{cl,^)R + ig^^-a^^a,^^ 

(2) 

where the potential V and the scalar-tensor coupling F 
are arbitrary real analytic functions. (The RGEs for this 
system have been studied earlier in [14] Although it is 
not necessary for some of the results of this paper, we 
shall assume that the potential has its minimum at the 
origin. Then, we can identify F{jS) = k = I/IGttG, G 
being Newton's constant, and V{Q) = 2kA, where A is 
the (dimension- two) cosmological constant. It will ap- 
pear that the behaviour of the couplings implicit in the 
functions V and F is sufficiently systematic that we are 
actually able to draw several conclusions involving an 
infinite number of couplings. For some purposes, how- 
ever, we shall restrict our attention to a five-parameter 
Ansatz, where V is at most quartic in </> and F is at most 
quadratic. 

Aside from establishing the existence of a nontrivial 
FP, the main question addressed in this paper will be 
the dimension of the critical surface. In practice, this 
is done by linearizing the flow around the FP. We define 
Vi ~ 9i^9i* the shift from the fixed point. The linearized 
flow around the fixed point is described by the equations 

dtVi = MijVj, (3) 

where My = Let P be the (generally complex) 

linear transformation that diagonalizes M: P~^MP = 
diag{ai, . . . , oat) (the columns of P are the eigenvectors 
of M). Defining / = P^^v, one finds dtfk = Oikfk, so 
fk{t) = e"''*. Transforming back to the original variables, 
the solution can be written gi{t) = +Re(Py/j(t)). It 
is easy to show that the eigenvalues ai are invariant under 
redefinitions of the couplings. 

The eigenvalues with negative real part (which for 
brevity we shall call the "negative eigenvalues" ) corre- 
spond to directions for which the RG flow approaches the 
FP in the UV. The corresponding parameter fi is called 
a relevant parameter. Those with positive real part (the 
"positive eigenvalues" ) correspond to directions for which 
the RG flow moves away from the FP in the UV. The cor- 
responding parameter /, is called an irrelevant parame- 
ter. The parameters corresponding to purely imaginary 
eigenvalues are called marginal. In the linearized theory, 
in order to approach the FP in the UV it is therefore nec- 
essary to stay on the hyperplane spanned by the eigen- 
vectors with negative eigenvalues. This hyperplane is the 
tangent space to the critical surface at the FP. Therefore, 
the dimension of the critical surface is equal to the num- 
ber of negative eigenvalues. On general grounds, one ex- 
pects the critical surface to be flnite dimensional. In this 
way all but a finite number of couplings would be fixed 
and the theory would be as predictive as a perturbatively 
renormalizable theory. 

We now give a brief summary of the results of this 
paper. First of all, a nontrivial FP still exists with the 
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Ansatz (2). The purely gravitational couplings (cosmo- 
logical constant and Newton's constant, which appear as 
the (/)^-independent terms in the functions V and F) have 
the same values as in [13] and all the other couplings are 
equal to zero. In a sense this is therefore "the same" 
fixed point that was considered in [8] and in [9, 15]. It 
can also be regarded as a generalization of the Gaussian 
FP of the pure scalar theory in flat space. We therefore 
call it the Gaussian-Matter FP (GMFP). We have per- 
formed a systematic search for other FP's within a five- 
parameter truncation of the action, where V and F are 
polynomials containing at most terms of order 0^ and 
0^ respectively. Thus, in addition to the cosmological 
and Newton's constant, we consider a scalar mass term, 
a quartic sclf-intcraction and a nonminimal coupling of 
the scalar field to the scalar curvature. Detailed numer- 
ical analyses have convinced us that there arc no FP's 
with nonzero scalar mass and couplings, for values of the 
cosmological and Newton constant close to the ones of 
the pure-gravity FP. 

Comparing to the results of the pure scalar theory, 
the main effect of the coupling to gravity is to change 
the exponents a^. It turns out that of the two canoni- 
cally marginal couplings, the (f>* coupling becomes irrele- 
vant while the (j>^R coupling becomes relevant. The other 
couplings preserve the character that is implied by their 
canonical dimension; the critical dimension would thus 
be equal to four. 

We then look at the effect of other matter fields on 
the FP. The results of this investigation generalize those 
already reported in [13]. The behaviour of the beta func- 
tions is determined by two parameters that depend on 
the number of fields, and the existence of the FP de- 
pends on the values of these parameters. In this way the 
existence of the FP yields constraints on the type and 
number of matter fields. These constraints appear to be 
satisfied by popular unified models. The existence region 
is subdivided into subregions with varying numbers of at- 
tractive directions. In the region that we have explored, 
comprising large numbers of matter fields, the dimen- 
sion of the critical surface is always finite. In particular, 
there are regions in which the attractive directions corre- 
spond to nontrivial polynomial potentials of degree four 
or higher. This yields a neat solution of a long-standing 
puzzle. In a pure scalar theory in flat space, the Gaussian 
FP is IR attractive (all couplings are irrelevant). As a 
consequence, when one takes the continuum limit at the 
Gaussian FP, the renormalized theory is free. This re- 
sult is not an obstacle in the context of an effective field 
theory, but it has to be somehow circumvented if scalar 
fields have to appear in a fundamental field theory. The 
coupling to gravity is a natural context for a solution of 
this issue. Our results imply that there exist theories 
of gravity coupled to matter such that the renormalized 
scalar potential has finitely many nonzero couplings in 
the continuum limit. This seems to indicate that the 
interaction with gravity (and, indirectly, with the other 
matter fields) solves the problem of the triviality of the 



scalar theory. 

This paper is organized as follows. In Section II, by 
way of introduction, we derive the ERGE and we use it 
to prove some well-known results on the Gaussian FP 
of a pure scalar theory in flat space. In Section III we 
consider the modifications of the beta functions due to 
gravity, we prove the existence of the Gaussian-Matter 
FP (GMFP) and we discuss the (negative) results of the 
numerical search for other FP's. Section IV is devoted to 
the properties of the GMFP. We analyze the linearized 
flow around the GMFP and show that the coupling to 
gravity affects the dimensions of the couplings, shifting 
them relative to the canonical values. In Section V we 
consider the effect of other massless, minimally coupled 
matter fields on the GMFP. In Section VI we will consider 
in some detail the dependence of our results on the shape 
of the cutoff function and on the gauge-fixing parameter. 
Finally in Section VII we make some concluding remarks. 

All results arc derived in the case of Euclidean signa- 
ture, in four dimensions. Since the expressions of the beta 
functions are extremely lengthy, in deriving our results 
we have made extensive use of algebraic manipulation 
software. 



II. THE GAUSSIAN FP IN PURE SCALAR 
THEORY 

We begin by considering the case of a single scalar field 
without gravity and a generic even potential 



1/(0) = £ A2 



i2n 



(4) 



In this section we assume that V{Q) = Aq = 2Ak = 0; the 
first nonzero term is the mass A2 = ^fn^ , while A4 is the 
usual quartic coupling. The couplings A2n have dimen- 
sion mass'^^^~^^\ so the usual power-counting arguments 
tell us that the terms in (4) with n > 2 are perturbatively 
nonrenormalizable, while the term rt = 2 is marginal. We 
will now rederive this result within the formalism of the 
ERGEs. This will set the stage for further developments 
in later sections. 

To derive the ERGE, one begins by modifying the clas- 
sical propagator by adding to the action a term quadratic 
in the fields which in momentum space can be written 
ASkicj}) = y dS(P{-q)Rk{z)c^{q), where z = q\ The 
effect of this term must be to suppress the propagation 
of field modes with momenta smaller than k, while leav- 
ing the modes with momenta larger than k unaffected. 
This is the case if the smooth cutoff function Rk is cho- 
sen to tend to zero for z k^ and to a constant for 
z 0. For numerical work, in this paper we will work 
with cutoffs of the form 



Rk{z) 



1 



-2az/fc2 



(5) 
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with a a free parameter. We then define a scale- 
dependent generating functional of connected Green 
functions 

Wk[J] = -In I (i?0)cxp(-(5 + AS, + J J<i>)) 



such that ^ 



= (0) and a scale-dependent effec- 
tk[(t)ci] - A5fe[0c;], where Vk[cj)ci\ = 



,7=0 

tivc action rfe[0c/] 

[J] — / J4>ci is obtained from Wk [J] by the usual 
Legendre-transform procedure. The scale-dependent ef- 
fective action tends to the bare action S when k tends 
to the UV cutoff, and to the ordinary effective action for 
A: ^ 0. We have 



dtWk 



dt{ASk) = ixr 



Rk 



(6) 



where the trace is over all Fourier modes (and internal 
indices, if there were any). Then, 



=dtWk[J] ~ dtASk 



2Tr((00) 



- ir- — —Of Hi.. 
2 5J5J * 



Applying the standard identity 

5^Wu ( 5^fk 



6J6J \ S(f>ciS 
one then obtains the ERGE [16] 
1^ fS^T, 



(7) 



(8) 



dtTk = -Tr 



+ Rk] dtRk. 



(9) 



In the previous formula and in the following we shall drop 
the subscript in (f>ci] this should not cause any confusion. 
The ERGE describes the flow of the functional Ffc with 
the scale k. In order to extract beta functions, one has 
to resort to approximations. A common procedure is to 



make an Ansatz about the form of and to insert it into 
the ERGE. Of course the beta functions obtained in this 
way are no longer exact: one loses all information about 
the dependence of the beta functions on the parameters 
that have been left out of the Ansatz. Nevertheless, the 
results do contain information that is not accessible in 
perturbation theory and they have been shown to yield 
numerically accurate values in many circumstances [4, 
17]. We now apply this procedure to the scalar theory. 

Introducing in equation (9) the truncation 
rfc(<^) =/d42; [-1^920 + y(02)j^ ^liere V is a k- 

dcpcndcnt potential, gives 



(10) 



where a prime denotes the derivative with respect to (fi^ , 
the trace can be understood as an integration over mo- 
menta. It can be reexpressed as: 



dtPk 



(11) 



327r2^^ \Pk + V + 4:(t>^V" , 
where Vol ~ J d'^x denotes the volume of spacetimc and 
1 



Qn[f] = 



Jo 



+ 00 



dzz^-'f{z). 



(12) 



The coupling constants can be extracted from the poten- 
tial by 



A2 



1 d'^V 



nl 9(02)n 



<t>=0 



(13) 



In order to look for a fixed point one has to define dimen- 
sionless couphngs A2„ = fc^("~^^2n- The corresponding 
beta functions are given by 



^,2(n-2) ^ Qn 

a.A2„ = 2(n-2)A2„ + ^-^ftr. 



4=0 



(14) 



Explicitly, the first few beta functions are given by 



9tA2 
dtX4 



dt^e — SAg 



" 32^^^ 
-3OA6Q2 
1 



dtp 



327r2 



dtXs =4A8 



327r2 
1 

32^ 
900A^Q2 



7' + 2A2 
dtV 

dtV 



144A^Q2 



-56A8Q2 

- 9OA10Q2 

dtV 
iV + 2A2)3 



dtV 
{V + 2A2) 



{V + 2A2)3 
72OA4A6Q2 



dtV 



1344A4A8Q; 



(V + 2X2 
dtV 



3 , - 1728A3Q2 



{V + 2X2Y 



(V + 2X2Y 



86A0XeXlQ2 



dtV 
{V + 2X2Y 



20736AIQ2 



dtV 



[V + 2A2)5 



(15a) 
(15b) 
(15c) 



(15d) 
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where Qn[f] ~ is a dimensionless integral, 

TZ = Rk is a dimensionless cutoff and and P = k~^Pk 
is a dimensionless modified propagator. 

Tliis theory admits a well-known Gaussian FP: if we 
set A2 = 0, equation (15a) implies A4 = 0, equation (15b) 
then implies Ae = 0, and so on: recursively all couplings 
are found to be zero. This is not the only solution of the 
coupled system. One can fix an arbitrary value of A2 and 
the equations then recursively determine all the other 
couplings [18]. However, when A2 7^ these potentials 
become singular at a finite value of (j) and therefore are 
not considered to be physically acceptable [19]. In what 
follows we will restrict our attention to the Gaussian FP. 

We now study the critical surface in the neighborhood 
of the FP using the linearized RG equation (3). Let 
P2n = dtX2n and let Mij = ff^. It appears from equa- 
tions (15) (as well as from dimensional and diagrammatic 
considerations) that the 2n-th beta function is a poly- 
nomial in the couplings A4, . . . , A2n-i-2, linear in A2,i+2- 
Therefore the elements of the matrix with j > i + 1 
are zero. On the other hand, since all A2n are zero at the 
FP, when the derivatives are evaluated at the FP only 
the terms linear in the couplings remain. These are ex- 
actly the terms on the diagonal, which are equal to the 
canonical dimensions of the couplings, and the terms on 
the second diagonal [i ^ j + 1), which are equal to 



M, 



i i+l 



df3: 



•2n 



8X2 



{2n+ l)(n+ l)c 



(16) 



where c = ^j^Q2 (^f ) All other terms are zero. Nu- 
merically, the integral Q2 (%?) is equal to 0.924 for 



Therefore, the matrix M has the following form: 



V 



-2 6c 

15c 

2 28c 

4 



(17) 



The eigenvalue problem for this infinite matrix yields the 
recursion relation 



A 



2n+2 



2(71 - 2) -yU 

(2n-hl)(n+l) 



-A 



2n 



(18) 



where /i is the eigenvalue. This relation can have two 
types of solutions. If we assume that the potential is 
a finite polynomial of order K, Eq.(18) implies that 
fi ~ 2{K — 2). These eigenvalues are just the diagonal 
elements of the matrix (17). The corresponding eigen- 
vectors are the columns of the following matrix P: 



1 




V... 



-0.0175512 3.84804 x 10 



-3 



0.999846 





-0.0438425 
0.999038 




1.04825 X 10"^ 
1.79148 X 10-3 
-0.0816446 
0.99666 



(19) 



The eigenvalues are equal to the canonical dimensions 
of the couplings, so that the relevant, irrelevant and 
marginal couplings correspond exactly to the couplings 
that are superrenormalizable, nonrenormalizable and 
renormalizable in the perturbative sense. 

These polynomial potentials suffer from the well- 
known problem of triviality. Consider the scalar theory 
regularized with a UV cutoff Aijv and the IR cutoff k. 
Keeping k fixed and letting Ajjv ^ 00 (the continuum 
limit) has the same effect as keeping A^y fixed and let- 
ting k 0. An irrelevant coupling tends to zero for 
A: — *■ 0, and therefore, for any fixed k it will tend to zero 
in the continuum limit. This will be the case for all A2i 
with z > 2, so the theory is non-interacting in the con- 
tinuum limit. (Our analysis only says that the couplings 
from Ag upwards have to be zero; the hard part is to 
prove that also the marginal coupling A4 tends to zero. 
For this, one has to go beyond the linearized analysis 
[20].) 

There is also another type of eigenvectors, correspond- 
ing to nonpolynomial potentials, that avoids the prob- 
lem of triviality. If we do not assume that A2_k-+2 = 
for some K, the recursion relation (18) can be solved 
for the A2n in terms of the free parameters A2 and /i, 
yielding a potential that can be written as a Kummer 
function [18]. There are (negative) values of ^ for which 
the potential has all the physically desirable properties 
(positivity at 00, symmetry breaking). They are there- 
fore nontrivial asymptotically free scalar theories. How- 
ever, there are infinitely many attractive directions and 
therefore these theories do not satisfy the conditions for 
asymptotic safety. 

This concludes our brief review of the ERGE for a 
scalar field theory. 



III. THE COUPLED SYSTEM 

We now consider the coupling of the scalar theory to 
gravity, using the Ansatz (2) for the running effective 
action. This will obviously change the beta functions 
of the scalar potential; in addition we will have to take 
into account also the beta functions of the gravitational 
couplings. These are given by the Taylor expansion co- 
efficients of the function F{(p'^) of equation (2), which we 
write as follows: 



(20) 



The first term in the expansion can be identified with the 
(inverse) Newton constant: £,0 ^ k ~ l/(167rG) while the 
second term is the well-known scalar tensor interaction 
term (p^R with dimensionless coefficient ^2 = The 
running couplings are given by 

~ 1 d''F 



6=0 



(21) 
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As 
by 



before, we define dimcnsionless couplings ^2n = 
"^•'^2n- The corresponding beta functions are given 



dt^2n = 2(n - 1)6 



1 



d 



Tl+l 



Vol n\ dRd{<PY 



dtVk 



p=0' 

(22) 

We now have to insert this Ansatz into the appropriate 
ERGE. The derivation of equation (9) in the previous 
section was quite general and therefore the ERGE for 
gravity coupled to a scalar field has again the same form, 
except for two generalizations: first, the field (p is to be 
reinterpreted as a matrix consisting of the components 
of the metric and a scalar field; second, since gravity is 
a gauge theory, one has to take into account the effect of 
gauge fixing and ghost terms. 

Here we mention some points that are necessary to un- 
derstand the results; we refer to [7] and [9] for details. In 
deriving the ERGE, one encounters the quantum metric 
(to be integrated out in the funtional integral), say 7^,^, 
which can be decomposed into the sum of an arbitrary 
background metric and a quantum fluctuation 
The background metric is used in the gauge fixing terms 
(23) below and also in the cutoff terms ASk, which have 
to be quadratic in /i^jy. In the Legendre transformation 
one encounters also the classical metric g^^, which is the 
canonically conjugate variable of the source associated 
to the quantum metric. Thus, in general, the action Tk 
will depend both on g and g. On the other hand the 
Ansatz (2) only depends on one metric. In order to de- 
rive the beta functions for the couplings in (2) we proceed 
as follows. In the r.h.s. of the ERGE, one first takes the 
functional derivatives w.r.t. the classical field g, then one 
sets the background metric equal to the classical one, i.e. 
9iii> = ff/^fi so that many contributions disappear. From 
these equations one can read off the beta functions of the 
couplings. 

The gauge-fixing action is chosen as 



(23) 



F„ 



-up 



^ugpp. 



/3 + 1 



so that the corresponding ghost action will be 



j2 fj,u 



(24) 

In principle, a and /3 are running parameters in the ef- 
fective action, so one should take into account their beta 
functions, too. However, as will be discussed in Section 
VI, there are arguments to the effect that a = at the 
FP. Therefore, unless otherwise stated, we will always 
work in the gauge a = and /3 = 1. 

The kinetic term of the gravitons is obtained by lin- 
earizing the action around a dc Sitter metric with scalar 
curvature R and a constant scalar backgound (j). Using 
the method of [7], the r.h.s. of (9) can be written as a 
sum of several terms, corresponding to the spin 2, 1 and 
components of the fields, and has to be completed by 
adding the ghost contributions. 

The spin-2 component of the metric has an inverse 
propagator 



1 



(25) 



where now z = — V^V^. The spin-1 component of the 
metric has an inverse propagator 



iFi^^)L+^-^R)-V{cb% (26) 



where a is the gauge-fixing parameter. The two spin-0 
components of the metric mix with the scalar field; the 
resulting inverse propagator is given by the matrix 



3 !3-a 
16 



R 
3 



2^ R\ , M,' Z + 2V' 



16 



iF'm 

iF'in 



^■A P, 



(z-|) + 0y' 



-|F'( 



(.-f)+0r 



^V" - R{2F' 



'■F") J 



(27) 



The two factors under trace in the r.h.s. of (9) are ob- 
tained from these expressions as follows. The modified 
(cutoff) propagators are given by the inverses of the ex- 
pressions in (25,26,27), with z replaced by Pk{z). The 
function Rk for each spin component is given by the dif- 
ference of the cutoff propagator and the original propa- 



gator. In the case of spin 2 and spin 1, this is just the 
function Rk{z) defined in (5), whereas for the spin com- 
ponents it is a 3x3 matrix (the difference of (27) with z 
replaced by Pk{z) and (27)). 

Since in (25,26,27) the momentum variable z always 
appears multiplied by the function _F(02), also the ma- 
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trices Rk appearing in (9) contain F[4P'). When inserted 
in the r.h.s. of (22), besides the explicit dependence of 
Pk{z) on fc, one has to take into account the dependence 
on k of all coupling constants that are present in F[(j)'^) 
and its derivatives (this is related to the "renormalization 
group improvement" that turns the one-loop RG into an 
exact equation). This generates terms proportional to 
the beta functions in the r.h.s. of the equations, so that 
the ERGE does not immediately yield expressions for the 
beta functions but rather linear equations for the beta 



J 



functions [24]. 

The beta functions themselves are then obtained by 
inverting the matrix of coefhcients, and this introduces 
further nonlincarities into the system. We will not write 
the expressions for the beta functions themselves but 
only the linear equations that determine the beta func- 
tions. We will order the couplings in order of decreas- 
ing mass dimension (before dividing by powers of k): 
Ao, A2, 6, A4, ^4 . . .. These are the first five equations: 



9tAo 



1 



327r2 
1 

3847r2 



dtViMiV + 8A2) + P(3P + 4A2)eo) 



dtV{'Xo{3V + IOA2) + V{11V + 26A2)eo) 



4o 



v{Xo - ioV) 



(28) 



+ Q2 
3847r2 ^0 



dtV{6Xo£,oV{20X2V + 20X1 + ^^(5 - 86)) + 3A2(-20A2-P - 20A2 + 7^2(8^^ _ 5))) 

7'2(7' + 2A2)2(67'-Ao)2 
'dtV{^ir^{-22QX2V - 220A1 +7^2(246 - 55)))" 



7'2(7' + 2A2)2(67'-Ao)2 



Qi 



7^(5Ao + 36P) 



Vi^oV - Ao 



5Q2 



7^(-3A2 + 6X0C0V + bilV 



2-D2M 



v^{ioV-x,Y 



(29) 



dtX2 



167r2 



Q2 



167r2^2-Q2 

5*6 r. 



dtV{2XlXi + 6^'(2A46 - A2(l + 26)") + Ao((l + 26)(-2A2 + 4A26^ + 6^') - 4A46^)) 

(7' + 2A2)2(67'-Ao)2 

7^(-A2(-4A2e2 + 8Ao667' + eo'^'(3 + 26)) + 6^(2A2 - 4Ao67' + (5 + 66)Co'^')) 



167r2^o 



7'(7' + 2A2)(67'-Ao)2 
n{-2Xo{V + 2A2) + 6^(5P - 2A2 + 1267'))" 



V{V + 2X2){ioV-Xo) 



(30) 



9*6 



487r2 



Qi 



487r2 



Q2 



(7' + 2A2)2(eo7'-Ao)2 

• {3A2A4 - Ao(4A2P(l - 36)6 + A^(3 + 166) + 7'(6A4Ao + V{1 - 36)6)) 
+ 6(10A^7' + lOA^ + 3A467'' + X2V^{1 - 66 - ^il))] 
dtV 



V^{V + 2X2f{^oV-Xo 



+ 



9*6 

3847r2e2 



Qi 



{-18Ai]7'2(-4A46 + {V + 2A2)6) - &Xl{-iXlV + 2A^ -|- XlV^{\ + 46 + 106^)) 
+ 6A2 (967^3 (-4A46 + ^iV) + A27'3(46 + 196' + 246' + 1866)) 
+ Ao6^(36A^ + SXlViQ + 76) - 12A27'3(6 + 186' + 186' + 966)) 
+ Ao67''(3A2(5 + 366 + 206') + 7'(2166A46 + 7^(106 + 216' + 366' - 5466))) 
- eo^'(104A^ + 3A2p2(23 _ 12^2) - 6A^7'(-25 + 46)) 
+ 26'A27'5(-5 + 246 + 516' + 366' + 1866) 

+ 36'7''(-24A466 - 76'P - 126'^^ + 666^) 

7^ 



P(P + 2A2)(6^-Ao)2 

{mXlilV + i2V{bXl - 10Ao67' + + 86)) 
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- 2X2i-5\lC2 + lOAoCo^^ + Co^'(-4 + 276))} 

n 



{lbXl{V + 2A2)'6 - 3AgCo(T' + 2A2)(8A^ + 156^' + 22A26^) 
+ ^o^^(-416A^ + A\lV{~92 + 616)) 

+ S.lV\£,2V{2?> - 1686 - 288C2) + 4A2(-20 + 796 + 60^^)) 
+ X^ilViUAXl ~ 20A26^'(-5 + 126)) 

+ AoCo'P^(4A^(18 + 376) + 6^^(85 + 1686 + 288^^))} 



5*6 



3847r26 

5*6 

3847r26 



Qi 

Q2 



7^ 



7'(7' + 2A2)(67'-Ao) 

n 



{~hXo{r + 2A2) + iioVi-V - I8A2 + 1667^)} 



p2(p + 2A2)2(67'-Ao)2 

{l^XliV + 2A2)' - 6Ao6^(12A^ + 8A27'(2 - 56) 



^2(5 + 286 + 966')) 



5* A. 



167r2 



Q2 



+ ^^7'2(-148A^ - 4A2P(31 + 606) + P'(-25 + 1686 + 576^2))} 
dtV 



(31) 



(7' + 2A2)3(67'-Ao)3 

• {Ag(-24A2 + h{V + 2A2)A6) + A§(-24A36 + 4Ai(8A4(l + 26) + V{-1 
+ A27'(72Al6 + 7'(-15A66 + UV{1 + 86))) 
+ 2A2A27'((-4A4(1 + 86 + 126') + 3(-5A66 + V{1 + 86)6))) 
+ 6(40A5 - 407'A^(-1 + 26) - 247'A^(-A46 + 567^)) 



166)6)) 



6'^'(24A26 - 5A667' + A47'(l 

f2 



86 + 246')) 

+ 6A2^'(-10A66' - 2A46(5 + 406 + 246') - ^(66' + 526' + 486^ 
+ Ao(20A4(1 + 46) - 2A3(-12A46 + V{1 + 246 + 1246' + 1266)))+ 
+ 2AoA27'(2A46(17 + 166) + 3P(6 + 46' + 446' + 266 + 16666)) 
+ AoA27''(-30A66' - 18A46(1 + 86 + 86') + ^(286' - 96^1 + 1466 + 6(1 + 6466))) 



866 + 6(1-1666))) 



5*6 

167126^ 



Q2 



Ao7''(72A^6'+^A46(l 

n 



86 + 246') + 7'(-i5A66' - ^(6 + 46' + 246^ - 66 - 8666)))} 



P(P + 2A2)2(6P-Ao)3 

• {(60A^6'P - 240A36'67'' + Xli-^XUl + 46V(9A4 + 76^^) - 24Ao6'^(6'^ + A06))) 
+ A^(8Ag6(36'^ + Ao6) 
+ p2(_2Ai] + 6A26P)(6' 



eo7'(7''(3 + 126 4 

- 66) 

+ 7'3(6'^(3A46(1 + 86 + 246') + 56' + 126' 
+ 3Ao6'^'(-A46(l + 86 + 246') + ^(-26' + 



3806') 



12Ao(3A4 + 6^))) 



f 726" - 566 - 24666)) 
366 + 8666)) 



- 2A27'((4Ag - 12A267')(6' - 66) - 6Ao6'^(A46(l + 86) + ^^(-26' + 66 + 4666))) 



167126'^' 



-2A2^i^7''(6A46(l + 86) + ^(86' 

n 



1326' -266 + 6(3 + 2466)))} 



7'(7' + 2A2)2(67'-Ao)2 

• {120A^6'P - 2A^(4A§6 - 8Ao667' + 6'^'(3 + 1906)) 
+ A27'(-8Ag6 + 16Ao667' + 486'^(2A4 + 6^)) 
+ eo(^'(3 + 166 + 3966') - 48Ao(2A4 + 6^))) 
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V\~2\l^2 ~ eoV{2A\M\ + 66) + V{i2 + 18^2 + 144^2^ - 24^0^4))) 
4AoCo^'(6A46(l + 66) + V{i2 - 666))} 



167r26^' 



(7' + 2A2)2(67'-Ao)2 



{-2Ao(7' + 2A2) + £,oV{bV - 14A2 + 246-P)} 



(32) 



This system of linear equations can be solved for /?2j = 
dt\2i and = 9t6i- Note that we have not made any 
truncation on the functions V and F: no couplings have 
been assumed to be zero. 

We do not exhibit the equations for the higher cou- 
plings. By means of algebraic manipulation software 
we have calculated the beta functions up to Ag and 6- 
The general pattern is however already clear from the 
equations shown here. The functions Qi and Q2 al- 
ways contain denominators involving only the couplings 
Ao, A2 and 6- Aside from these couplings appearing 
inside Qi and Q2, the equation for /3^j is a polyno- 
mial in the couplings Aq, A2, . . . , A2j+2 and 6j 6, ■ • ■ , 6i, 
with the functions Qi and Q2 as coefficients, while 
the equation for is a polynomial in the couplings 
Ao, A2, . . . , A2i+2 and 6: 6, • ■ • , 6i+2, with the functions 
Qi and Q2 as coefficients. When these equations are 
solved, the beta function /3^j is written as a rational 
function of Ao, A2, . . . , A2i+2 and 6: ?2, ■ . . , 6i: while the 
beta function [i^i is written as a rational function of 

•^0 1 -^2 : ■ • ■ , -^24+2 and 6 7 6 7 • ■ • J ^21-1-2 • 

It is clear that the system /Sji = 0, — admits a 
FP for which all couplings A2i and 6i vanish for i > 0, 
while for i = Aq = 2kA and 6 = have the same 
values that they would have in the presence of a single 
free scalar field, as discussed in [13] (these values arc 
numerically very close to those of pure gravity, discussed 
in [9]), namely, for a = 2, 



Ao, = 0.0080022, 
6* = 0.023500. 



(33) 



To compare with the results of [9] , we define the dimen- 
sionless variables A = A/fc^ = Ao/26 and g = Gk^ = 
1/I67r6- At the GMFP 



A, — 



Ao* 



= 0.1703, 



26* 

1 

5, = TT— = 0.8466. 



(34) 



16^6 



0* 



These values differ from those in eq. (5.25) of [9] on two 
accounts: they are calculated for different values of the 
cutoff parameter a, and here the FP is shifted due to 
the presence of the scalar field. When these factors are 
taken into account there is perfect agreement. (See figure 
(5.b) of [9] for the dependence of results on s = 2a and 
compare with figure 8.) 

This FP can be viewed alternatively as the FP of pure 
gravity, slightly shifted due to the presence of a free, 



massless, minimally coupled scalar, or as the Gaussian 
FP of the pure scalar theory, generalized to include grav- 
itational interactions. It is remarkable that matter re- 
mains "non-self-interacting" at this FP, and that the only 
nonzero couplings are those that affect only the gravita- 
tional degrees of freedom. (This goes some way towards 
justifying the assumption in [13] that matter fields are 
non-self-interacting.) For want of a better terminology, 
we shall refer to this FP as the Gaussian-Matter FP, of 
GMFP for short. 

The issue arises whether the coupled system of equa- 
tions admits other nontrivial FP's. The complexity of 
the equations has prevented us from deriving definite re- 
sults on this issue. We have looked for other FP's using 
numerical methods in a five-parameter truncation of the 
theory containing the couplings A2n for n = 0,1,2 and 
the couplings 611 for n = 0, 1. Our method consists in 
considering a grid in the space of all parameters, evaluat- 
ing the beta functions by numerical integration at a point 
and then at all neighboring points. If all beta functions 
change sign simultaneously when going from a point to a 
neighbor, then generically there will be a FP somewhere 
near the link between the two points. The region is then 
examined with a finer grid until the position of the FP is 
located with sufficient accuracy. We have started oS with 
the 2x2 grid given by Aq and 67 confirming the results 
of [9]; we have then added one by one the other vari- 
ables involved in the five-parameter truncation, getting 
increasingly complicated systems of equations. 

Because of the complexity of the beta functions, 
the numerical evaluation takes considerable time. The 
largest range we have explored is a 5 x 5 grid with the di- 
mcnsionless cosmological constant A ranging from 0.010 
to 0.045 in steps of 0.005; the dimensionless Newton con- 
stant g ranging from 0.01 to 0.06 in steps of 0.01; the 
dimensionless scalar mass 2A2 ranging from -1 to 1 in 
steps of 0.2; the dimensionless quartic scalar coupling A4 
ranging from -15 to 5 in steps of 1; the dimensionless 
scalar-tensor coupling ^ = 6 ranging from -5 to 5 in 
steps of 1. This makes a lattice with more than 120,000 
points. Many other attempts have been tried with finer 
lattices and/or fewer parameters. We did find some non- 
trivial solutions when considering less than 5 parameters, 
but none of them survived the addition of an extra pa- 
rameter, so that we had to conclude that they were all 
spurious FPs due to the truncation. As a further check 
we have also resorted to series expansions around the 
GMFP or one of the spurious FP mentioned above. All 
the results obtained in this way are perfectly consistent 
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with the other calculations. The outcome of all these 
efforts is that no FP other than the GMFP was found. 

This is of course not a proof that an FP docs not exist 
in this range of couplings. For example, a beta function 
may change sign twice on a link, once from positive to 
negative and once from negative to positive, and if the 
distance between the zeroes is smaller than the size of 
the step, it may well escape detection by our methods. 
Nevertheless, after this numerical work, we consider it 
quite unlikely that another nontrivial FP exists in the 
range of values for the couplings that we have considered. 

This result is corroborated by the following observa- 
tion. If one does not truncate the functions V and F to 
polynomials, as in the pure scalar case the structure of 
the beta functions seems to allow for a recursive solution 
depending on two free parameters. If we fix arbitrary 
values for Aq and ^o- from equations (28,29) one derives 
A2 and ^2; substituting them into equations (30,31), one 
can solve for A4 and ^4 and so on. This will determine 
the functions V ~ and F ~ F^, up to two arbitrary 
parameters. It will be interesting to analyze this in de- 
tail and to see whether the resulting functions and 
F^ are regular or still present the problems discussed in 
[19]. In any case, it seems highly unlikely that the solu- 
tions will be polynomial. This point of view also sheds 
a different light on the FP found in [8, 9]. The values 
of Ao* and ^0* at the GMFP are the only ones for which 
A2* = and ^2* = 0, and as a consequence all the higher 
couplings turn out to be zero, in accordance with the 
truncation made there. We shall not pursue this issue 
anymore here. In the rest of this paper we shall restrict 
our attention to the GMFP, which is a special member 
of this family of solutions, and is definitely a physically 
acceptable solution. 



We begin by calculating the matrix My. Again, we or- 
der the couplings in order of decreasing mass dimension: 
•^0 , $0 , A2 , , A4 , ^4 , Ae , . . . . As in the pure scalar theory, 
due to the functional dependences of /J^^ and on the 
couplings, an infinite triangle above the diagonal is zero. 
Furthermore, due to the fact that only the "purely grav- 
itational" couplings Ao and ^0 are nonzero at the GMFP, 
an infinite triangle below the diagonal is zero. The struc- 
ture of the matrix M is therefore remarkably simple: 



/ Moo 


M02 













M24 











M44 


M46 











M66 


V ••• 









(35) 



where each one of the nonzero entries is a 2 x 2 matrix 
of the form 



0X(2i) 



(2j) 



(36) 



9e 



For the calculation of the dimension of the critical sur- 
face we need to count the number of negative eigenvalues 
of the matrix M. The eigenvalue problem for the ma- 
trix Al could be turned into recursion relations for A2i 
and ^2i, as for the pure scalar theory. However, if we 
restrict ourselves to solutions where V and F are poly- 
nomials, given the almost-block-diagonal structure of M, 
the eigenvalues of M are just the eigenvalues of the diag- 
onal blocks Ma . Explicitly, the diagonal blocks have the 
following form: 



IV. LINEARIZED FLOW AROUND THE GMFP 

Having established the existence of the GMFP in the 
truncation defined by the action (2), we have to study 
its properties, in particular to find the dimension of the 
critical surface. 



M, 



2{i - 2) 
2{i-l) 



6Rhx 6Mx^ 



(37) 



where the first term contains the canonical dimensions of 
the couplings and the second term, which contains the 
quantum corrections, has the following form: 



<5Maa = 



J2_ 
167r2 



3Qi 



dtV 



(Ac - V^oY 



1 

A2' 





n 




-3Q2 







5Q2 



-8Q2 



7^(-2Ao + h£,^V) 
n{-2Xo + 5^oV 



dtV 



Vi-Xo + ^oV) 



Q 



(Ao - CoVr 

n 



(Ao - ^^Vf 



10^oQ2 
+ ioeoQ2 



7'2(Ao-^o7')2 

dtV 



i^oV - Ao)3 

n 



•A 



Qi 



Vi-Xo + ^oV) 



5Q2 



{^oV-Xo)\ 

vHXo-iovy- 



(38) 
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1 „ 



dtV 



(Ao - ^ovy- 

2 



327r2A2 



6M, 



7^ 

eAoCo • A • Q2 

2Q2 



^oV - Ao 
7^ 



(Ao - CoVr 



Qi 
Qi 



Ot7^(-3Ao + llgo7^) - 



5Q2 
-5Q2 



\n- 


2Ao 


+ 5^07')] 


■ (Qi 


IdtVi- 


-3Ao + 




[ n 


-Ao 




[ V{ 


-Ao + 








'dtP{5Xl 


- lOAoeoT' - 




)1 


-5Q 






^(Ao- 









5Q2 



dtvisxl-exo^ov + miv^) 



v^iXo-Covy 



3-D3M 



Qi 



- 2 • A • Q2 

9t7'(3A2-22CoAo7' + lieo'7'' 



dtV 



5Q: 



V^{-X,+^^Vf 

7^(-2Ao + 5go7^) 
^(-Ao + Co^) 

dtV{-^xl + 9Xl^oV - ssXoeoV + neoV^) 



4^0 (Qi 



(Ao - eoT')^ 
^(^(-SAo + lieo^) 



7'2(-Ao+eo7')3 



A 



7'(-Ao + Co7') 



(-Ao + eo7')3 

9t7'(3A2 - 6Ao6^ + 11^0^') 



-5Q2 

n 



(Ao - ^ovy 



n 



lOQ: 



(-Ao+^oP)^ 



(39) 



(40) 



Qi 



9t7'(-3Ao + lieo7') 



P(-Ao + eoP) 
7^(5Ao - lOAo^oT' - ie^V^ 



5Q2 



at7'(3A2-6Aoeo^ + lUo'^2) 



^(Ao - ^0^)2 
dtV{-iXn + ll^oT') 



7'(-Ao + ^o^) 
dtViiXl - 22AoCo7' + llCgp2) 



-5Q2 
5Q2 



V^{Xo-£.aVf 
7^(3A3 - 9XlioV + ITAo^o'^' + 5^o^') 



p2(-Ao + eoP)3 
at7'(3A2 - eAo^oT' + llQV^ 



•A 



where 



A = Qi 



7'(Ao - $0^)2 

7^(5Ao + 3eo7') 



-5Q 



7'(Ao - i^V) 



5Q2 



7'2(Ao - ioVf 
dtVi-SXl + QAg^oT^ ^ 33Aoggp2 + llCgT^^) 

^2(-Ao + ^0^)3 

7^(3A2-6Aoeo7'-5eo27'2)^ 



-p^Xo-^ovy 



3847r^$o 



The most remarkable property of these quantum cor- 
rections is that they are independent of z, so that the 
eigenvalues of M2i2i simply grow by 2 whenever i is in- 
creased by 1. For example, choosing the cutoff with 
a — 2, we have the following numerical results: 



I 

with eigenvalues — 0.1173 ± 3.1563i; 

_ / 5.1257 -2.5192 
AUi - 1^ 8.1295 -1.3604 

with eigenvalues 1.8826 ± 3.1563i; 



Mm 



1.1257 -2.5192 
8.1295 -5.3604 



which has eigenvalues —2.1173 ± 3.1563i; 

3.1257 -2.5192 \ 



M22 = 



^.1295 -3.3604 J 



(43) 



(44) 



7.1257 -2.5192 
8.1295 0.6396 



(41) 



(42) 



(45) 



(46) 



with eigenvalues 3.8826 ± 3.1563i, and so on. 

The off-diagonal blocks Ma^i in (35) do not affect 
the eigenvalues but determine the mixing between the 
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couplings. Numerically, we have 



02 



M24 = 



M46 = 



-0.005036 
0.002736 



-0.002264 
-0.007585 



-0.03021 -0.01359 

0.01642 -0.04551 

-0.07554 -0.03340 

0.04104 -0.1138 



(47) 



(48) 



(49) 



and so on. 

The first two (complex conjugate) eigenvectors have 
components 

/ -0.3486 ±0.3392i \ 
0.8737 



. (50) 







V 



They are a mixing of Aq and ^0; the corresponding 
(complex conjugate) eigenvalues have negative real part 
—2.1173 and therefore these are relevant couplings. 
The second and third eigenvectors have components 

/ (-16.59T 5.343i) x 10"^ 
(-2.970 ± 1.136i) X 10"^ 
-0.3485 ±0.3392i 
0.8736 







V 



(51) 



They are essentially a mixing of A2 and ^2, with 
small contributions from Aq and S,o', the corresponding 
(complex conjugate) eigenvalues have negative real part 
—0.1173 and therefore also these couplings are relevant. 
Since they lie very close to the plane spanned by the 
relevant coupling A2 (canonical dimension 2) and the 
marginal coupling ^2 (canonical dimension 0), we can say 
by a slight abuse of language that the quantum correc- 
tions change the dimension of A2 and ^2 making them 
both relevant. 

The fifth and sixth eigenvector have components 



(19.34 T2.834i) x 10"*^ 
(2.653 T2.310i) x 10"^ 
(-9.95lT3.203i) x lO^^ 
(-17.81 ±6.820i) X lO^^ 
-0.3485 ±0.3392i 
0.8727 





(52) 



They are essentially a mixing of A4 and ^4, with small 
contributions from Aq, ^Oi A2 and ^2; the corresponding 
(complex conjugate) eigenvalues have positive real part 
1.8826 and therefore these couplings are irrelevant. Since 
they lie very close to the plane spanned by the marginal 
coupling A4 (canonical dimension 0) and the irrelevant 
coupling ^4 (canonical dimension -2), we can say by a 
slight abuse of language that the quantum corrections 
change the dimension of A4 and ^4 making them both 
irrelevant. 

The pattern continues. The eigenvalues come in com- 
plex conjugate pairs, and are formed by mixing the cou- 
plings X2i and ^2i, with small contributions from the 
lower couplings. The eigenvalues also occur in complex 
conjugate pairs, and arc equal to (minus) the canonical 
dimensions of the couplings {2{i — 2) and 2(i — 1) re- 
spectively) plus a quantum correction. The correction is 
positive for the couplings of the series A2i and negative 
for those of the series ^2i, and the resulting dimension 
is always contained between those of the two main cou- 
plings that enter into the mix. 

All eigenvalues differ from the first two by multiples of 
2. In particular, all the eigenvalues from the fifth onward 
have positive real parts, so that the dimension of the 
critical surface is four. The naive expectation based on 
canonical dimensions would have been five (or three, if we 
don't count the two marginal couplings). The quantum 
corrections modify the dimension of the two marginal 
couplings A4 and ^2 so that ^2 (after mixing with A2) 
becomes relevant while A4 (after mixing with ^4) becomes 
irrelevant. 



V. EFFECT OF OTHER MATTER FIELDS 

In this section we assume that in addition to the gravi- 
ton and the scalar field discussed in the previous sections, 
there are ng- 1 new real scalar fields, nw Weyl fields, um 
Maxwell fields, n^s (Majorana) Rarita-Schwinger fields, 
all minimally coupled. We neglect all masses and inter- 
actions of these additional matter fields. The only inter- 
actions arc the ones discussed in the previous sections. 
This generalizes the results of [13] where only the cou- 
plings Aq and ^0 were taken into account. We also give 
some more details of the calculations. 

In the presence of these new fields, the equations 
(28,29) for the beta functions are modified by the ad- 
dition of the following terms: 



SdtXo 



Sdt^o = 



1 



327r2 
1 



{ns - 2nw + SriM - 4n_Rs) Q2 



dtV 



2 . {-2ns + ^nw - ^um) Qi 
dtV 



3847r 

+ {-6nw + QriM - 1671 i?,s) Q2 



dtV 



V 



(53) 



(54) 



Since the contribution of the new fields to the effective 
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action is independent of cj), they do not affect at all the 
beta functions of all couplings and for i > 1. 

The equations for the couplings \2i for i > 1 and for 
i>\ are automatically satisfied at the GMFP. Therefore 
the only equations that remain to solve are the ones for Aq 
and ^0- For the sake of comparison with [13] we will use 
the couplings A and g in place of Ag and ^o- The system 
of these two equations is the same as the one discussed in 
[13], and therefore the values of Aq and at the GMFP 
coincide with the ones calculated therein. 

We recall some of the calculations in [13]. For the 
purpose of finding the fixed points, one can use the fol- 
lowing trick. We observe that at a fixed point (?t^o/Co = 
—dtg/g = 2. Therefore, the equations for the fixed points 
are equivalent to another, simpler, set of equations which 
is obtained by replacing 9f^o/^o with 2 in the r.h.s.'s of 
(28,29). Then, the equation = can be replaced by 



g ■ c(A) - 2A = 0, 



(55) 



where c(A) is obtained by formally replacing G with 1 
and dtG with —2 in the expression for dtA/k^. When 
c(A) ^ 0, we can solve (55) for g and substitute the result 
into (3^ = 0. We shall denote 



h{X) - /3M A, 



2A 

4X) 



(56) 



so that the zeroes of h correspond to the FP's. 

The general behaviour of the function h is controlled 
by the values of two parameters A' and r, which in turn 
depend on the type and number of matter fields. The 
parameter A' is equal to A + cr, where A = n;, — n/ is the 




FIG. 1: Plot of the function h for ns = 0,nw = 0,nM = 
(Region I); the asymptotic value of /i for t — > — oo is -65.3. 



difference of the total numbers of bosonic and fermionic 
degrees of freedom {rif = 2nw + 4n and nb = ns + 
2nM + 2) and cr = 20 Q2 [^] /Q2 [^] is approximately 
equal to 3.64 (for a = 2). 

The value of the cosmological constant at the FP, A*, 
is zero on the hyperplane A' = 0. To see this, note that 
when c(A) = 0, equation (55) implies A=0. Therefore, if 
c(0) 7^ the only solution with A* = is the Gaussian 
FP, but if c(0) = we can have a GMFP with A, = 0. 
Explicitly, 



c(0) 



1 



47rfc4 



{nb ~ n/)Q2 



dtV 



-2OQ2 



(57) 



so that the condition for the existence of a non-Gaussian 
FP with zero cosmological constant is precisely 



A' = 



(58) 



Due to the irrationality of cr, there is in general no combi- 
nation of matter fields that satisfies this condition; how- 
ever, the hyperplane defined by (58) has an important 
physical significance: it separates the regions with posi- 
tive and negative A*, as shall become clear below. 

The function h tends to zero when A — > 
minze[o,oo] (^)/2. However, this point does not corre- 
spond to an FP: For this value of A the denominators 
in the functions Qi and Q2 appearing in the beta func- 
tions vanish and the beta functions themselves blow up. 
Moreover, it was shown in [6] that the Ward identities 
break down near this point. Consequently, only values of 
A strictly less than minj,g[o,oo] (^)/2 will be considered. 

The function h(X) always has a zero at the origin, cor- 
responding to the Gaussian FP. The derivative of h{X) 
at the origin is given by 



IGtt 



1 

A' 



(59) 



and therefore has the same sign as A'. When A' > 0, the 
function h(X) tends to —00 for A somewhere between 
and min^g[o.oo](^)/2 (namely where c(A) = 0). Conseqe- 
untly, there exists a non-Gaussian FP with positive A*. 
On the other hand, when A' < 0, h has no positive zeroes 
and the existence of the NGFP for negative A, hinges on 
the asymptotic behaviour of h for A — ^ — 00: it only ex- 
ists if h tends to a negative asymptote. The asymptotic 
behaviour of h is given by lim ^(A) = 1927r/T, where 

A — ^ — 00 

T =To + nsTs + nwTw + nMTM + njisTRS (60) 
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FIG. 2: Plot of the function h for ns = 100, nw = 40, um = 
(Region II). The asymptotic value of h for t — oo is 0.939. 
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FIG. 3: Plot of the function h for ns ~ 0,nw ~ 40,nA/ = 
(Region III). The asymptotic value of h for t —oo is -8.13. 



tm =Qi 



dtV 



r 



9Q 



dtV 



trs =16Q2 



dtV 



p2 

14.79. 



-6.52 



Depending on the sign of the two parameters A' and 
T, the space spanned by the variables ns, nw, nM and 
niis can be divided into four regions that we shall label 
as follows: 





T < 


r > 


A' 


< 


III 


IV 


A' 


> 


1 
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(The numerical values are given for a = 2). 



The behaviour of the function h is shown in figure (1) 
for pure gravity, which lies in region I. There are no ze- 
roes for negative A, since h grows monotonically from the 
asymptote at t — oo to zero in the origin (Gaussian 
FP). It then has a positive zero and tends to — oo. There 
is another apparent zero for A ~ 0.4, but it is not an ac- 
ceptable solution: it corresponds to the point where the 
denominators in the function Qi vanish. Thus, in region 
1 there is always a single FP with positive A*. 

Figure (2) shows the function h for a theory in region 
II. The behaviour for positive A is very similar to that in 
region 1, but the asymptote for i — ^ — oo is now positive, 
so that there exists a second FP for negative A* This 
FP can be seen to yield negative g^, and is therefore 
physically uninteresting. 

The behaviour of the function h in region 111 is shown 
in figure (3). The positive zero is the unphysical one, 
so there is a single attractive GMFP with negative A*, 
which turns out to have positive . 

Finally, the behaviour of the function h in region IV 
is shown in figure (4). It decreases monotonically from 
the positive asymptote t — > — oo to the Gaussian FP. For 
positive A it behaves like in region 111, having no zeroes 
except for the unphysical one. Thus, in region IV there 
is no non-Gaussian FP. Region IV is the white wedge 
in figures (5,6,7). One sees that it comes actually quite 
close to the origin; from this point of view the existence 
of the FP for pure gravity seems to be a lucky accident. 

The value of A* in regions I and II is always less than 
min2g[o,oo](^)/2, which is numerically equal to 0.402 (for 
a = 2) and therefore reasonably within the bounds of the 
heat-kernel approximation. On the other hand in region 
111 A* becomes quickly rather large in absolute value; in 
this regime on shell and therefore the heat-kernel 

approximation ceases to be valid on shell. In this region 
the results are only reliable close to the surface A' = 0. 

In order to determine the dimension of the critical sur- 
face we have calculated numerically the matrix M for 
many different combinations of fields. The results of 
such calculations are shown in figures (5,6,7) for the case 
nRs = and nM = 0, riM = 24 and nM = 45 respec- 
tively (these numbers are chosen to correspond to the 
gauge field content of popular GUT models). 

First of all, these numerical calculations exactly con- 
firm the shape of the existence region of the FP that was 
derived analytically above. The structure of the eigen- 
values is the same as in the pure gravity-|-scalar case, 
which was discussed in Section IV. The eigenvalues are 
given by the canonical dimensions plus a quantum cor- 
rection which depends on the type and number of matter 
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FIG. 4: Plot of the function h for ns = 40, nw = 40, um = 
(Region IV). The asymptotic value of h for t — ^ — oo is 8.72. 



fields but otherwise is the same for every pair (A2i,^2i)- 
For any given number of matter fields, the GMFP has a 
finite-dimensional critical surface. In region III, the crit- 
ical surface has mostly dimension 3 (three negative real 
eigenvalues), except for a narrow area close to the sepa- 
ratrix A' = 0, where its dimension is 2 or 4. In regions 
I and II the critical dimension varies considerably, being 
roughly linear in the number of fields (it grows with ns 
and decreases with nyy)- These calculations generalize 
the results of [13], where the FP could have at most 2 
attractive directions. 

It is interesting to compare these results with the anal- 
ysis of the pure scalar theory. Due to the fact that the 
coupling A4 is marginal in the pure scalar theory, the 
linearized analysis is not sufficient to determine its be- 
haviour. In the presence of gravity there is no zero eigen- 
value and therefore the linearized analysis is sufficient to 
determine the dimension of the critical surface. In re- 
gion III, the relevant directions correspond to potentials 
that are at most quadratic in (/). In region II, however, 
there can be a large number of negative eigenvalues, cor- 
responding nontrivial potentials V that are polynomial 
and asymptotically free, thus avoiding entirely the trivi- 
ality issue. These theories are also predictive since they 
have a finite number of negative eigenvalues [25]. They 
are therefore asymptotically safe. 



VI. CUTOFF AND GAUGE DEPENDENCE 

Physical results are independent of cutoff parameters 
in the exact theory, so the extent of parameter depen- 
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FIG. 5: Attractivity regions for njv/ = 0. The colours corre- 
spond to the number of attractive directions. 
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FIG. 6: Attractivity regions for um = 24. The colours corre- 
spond to the number of attractive directions. 
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FIG. 7: Attractivity regions for um = 45. The colours corre- 
spond to the number of attractive directions. 



dence that is observed in the truncated theory gives a 
quantitative measure of the errors. We have performed 
various tests on the parameter dependence of our results 
and it is reassuring for the reliability of the truncation 
that this dependence turns out to be reasonably mild. 

The dependence of Aq* and ^0* on gauge and cutoff pa- 
rameters was discussed in [9, 21]. Figure (8) summarizes 
the cutoff-dependence at the GMFP for gravity coupled 
to one scalar field in the gauge a = 0. The results we ob- 
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FIG. 8: a-dependence of Ao*, Co* and Ao*/Co* in the gauge 
a — 0. The values of Ao* and ^o, have been magnified by a 
factor 500 to display the three curves in the same range. 
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FIG. 9: a-dependence of the real part of the eigenvalues of 
the stability matrix (first 2x2 submatrix). 



tain are very close to those of [9] , since at the GMFP the 
only new contribution that we get for the values of Aq* 
and ifo* is that of the kinetic term of the scalar field. It is 
apparent that while Aq* and ^o* are quite sensitive to the 
cutoff parameter a, the ratio Xo*/£,o* is not. As noted in 
[7], this quantity is, up to numerical factors, (the inverse 
of) the on-shell action, a physically observable quantity, 
so it must be independent of the cutoff scheme. It is seen 
in figure (8) that its a-dependence is indeed pretty mild. 

The dependence of ± 162^, the eigenvalues of the 
stability matrix M, on the cutoff parameter a is shown 
in figures (9-11), for several values of a. We have cal- 
culated them in the range 1/5 < a < 50, but they are 
only reported for 1/2 < a < 20. Figure (9), giving the 
real parts of the eigenvalues of the submatrix Afgo i agrees 
with figure (9) of [9], up to the small corrections due to 
the presence of a scalar field. The figures relative to real 
parts of the remaining eigenvalues are simply shifted by 
the canonical dimension 2i. 

The first thing we can notice in figure (9) is the pres- 
ence a clear plateau with very weak (apparently logarith- 
mic) variation of the eigenvalues, for 1 ^ a < 20. Actu- 
ally, the results for 1/5 < a < 1/2 seem to indicate that 
there is a divergence as a ^ 0. This is due to the fact 
that in this limit the cutoff function tends to become a 
constant, so it affects also the propagation of modes with 
momenta larger than fc, and it does not work well as in 
IR cutoff. Clearly, larger values of a, of order unity, are 
preferred. This is in accordance with the generic features 
of the cutoff functions described in [5] . 

As far as the a dependence is concerned, we can notice 
that it is quite weak. For all possible values of a all curves 
are contained between the curves a = and a ~ 3, 
which differ by ~ 0.4. In order to better understand the 
dependence on a for different values of a, it is useful to 
plot the same results as a function of a (figure (10)), 
with a being a parameter that labels different curves (we 
shall restrict only to one plot of the real parts, the others 
can be obviously derived by shifting the graph by the 
canonical dimension of the operator involved) . 

Since the real parts of the eigenvalues of the matrix 
M22 are close to zero, this modest shift of the eigen- 
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FIG. 10: a-dependence of the real part of the eigenvalues of 
the stability matrix (second 2x2 submatrix). 



values due to the change of gauge parameter is enough 
to change their sign. For instance, in the gauge a = 1 
one would compute the dimension of the critical surface 
to be two. Looking at figure (10) one can most easily 
understand what the situation is like: for a = all cut- 
offs give a negative value of 9' , then as a increases they 
change sign, but for large values of a they become nega- 
tive again. Physical results such as the dimension of the 
critical surface cannot depend on the shape of the cutoff 
function, so this fact is certainly a shortcoming of our 
truncation. More work is needed to assess with greater 
confidence the dimension of the critical surface, but the 
considerations developed in [9], i. e. that a itself runs to 
in the UV regime, suggest that the a = value is the 
physicaly correct result. 

The same conclusions for the cutoff and gauge inde- 
pence can be drawn for the imaginary parts, as can be 
seen from figure (11); they turn out to take the same 
values (up to a sign) for all the eigenvalues. The effect 
of nonvanishing imaginary parts is that the RG spirals 
around the FP, but they are not important in the discus- 
sion of the attractivity of the FP. 

This discussion applies also to the higher couplings; 
their a- and a-dependence is given by curves that dif- 
fer from those in figures (9-11) by a constant shift by a 
multiple of 2. The only important point that remains to 
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FIG. 11: a-dependence of the imaginary part of the eigenval- 
ues of the stabihty matrix. 



some extent open, then, is the exact dimension of the UV 
critical surface, but nevertheless we can safely say that 
it is finite-dimensional. 

When other matter fields are present, the nature of 
the GMFP as a function of the number af matter fields 
is also a function of gauge and cutoff parameters. As 
already noted in [13], the constant a is independent of 
the gauge parameter and varies from 4.745 for a = 0.05 
to 2.765 for a = 20. This corresponds to a vertical shift of 
the separatrix A' = by at most 2 in figures (5-7). The 
parameter r is gauge independent for a ^ 0, but shows 
discontinuity at a = 0. The plane r = is shifted and 
also slightly rotated to the right as a grows. Thus, region 
III becomes larger as a grows. Recall however that only 
the part of this region close to the separatrix A' = is 
trustworthy. 



VII. CONCLUSIONS 

In this paper we have considered the application of 
the ERGEs to a coupled system of gravity and matter 
fields. The main aim of this work was to verify that the 
conditions for asymptotic safety continue to hold in the 
presence of interacting matter fields. To make the prob- 
lem manageable, we have first dealt with a single scalar 
field (j) with an arbitrary potential depending on 0^, to 
see how this inclusion could change the picture of pure 
gravity, then we have considered the effect of minimally 
coupled fields with different spins. Our results can be 
considered as a first step towards constructing a realis- 
tic theory of gravity and matter, but are also relevant to 
gravitational theories containing a dilaton. 

In the context of the Ansatz (2) we found that there 
exists a FP where only the cosmological constant and 
Newton's constant are nonzero. We called it the "Gaus- 
sian Matter" FP. A detailed numerical search within a 
five-parameter truncation of the effective action has failed 
to yield any other FP. This is actually what one would 
expect from our understanding of the scalar theory [18]. 

The GMFP may be viewed in two ways. On one hand, 
the scalar field can be regarded as a "perturbation" of the 



pure theory of gravity considered in [9] and the GMFP as 
an extension of the FP found in [8] . The addition of the 
scalar field has the effect of shifting slightly the values 
of A and already noticed in [13], while the attrac- 

tivity property is preserved. On the other hand, we can 
consider the effect of adding gravity to a scalar theory 
with a generic potential and regard the GMFP as an ex- 
tension of the Gaussian FP. The main effect is that the 
couplings that are present in the scalar potential V mix 
with those appearing in the function F, and the dimen- 
sions of the resulting couplings (which dictate the speed 
of the approach to the FP) is changed by a finite quan- 
tum correction. While at the Gaussian FP the quantum 
corrections vanish, so that the relevant couplings are, as 
usual, those with dimension less than four, the gravita- 
tional contributions bring about modifications even if the 
matter sector allows for a perturbative treatment. At the 
GMFP, the stability matrix has a block-diagonal form so 
that there is a strong ing between the parameters X2n and 
S.2m for n ~ m, whereas for n ^ m they are almost or 
completely decoupled. The eigenvalues, whose real part 
determines whether an operator is relevant or irrelevant, 
come in complex-conjugate pairs, and grow systemati- 
cally by a constant 2. For instance, the marginal operator 
of the pure-scalar cj)^ theory becomes now an irrelevant 
operator, and the dimension of the UV critical surface is 
calculated to be 4 for a generic analytic potential. These 
results hold in the gauge a = 0; they differ slightly for 
other values of the gauge-fixing parameter, but a ~ 
seems to be the physically correct value at the FP. The 
striking fact is that gravity gives calculable, finite contri- 
butions that change significantly the pure-scalar theory. 
This is one of the most important results of our paper. 

We have then considered the effect of adding other min- 
imally coupled massless matter fields. For the existence 
of the GMFP, we obtain the same bounds presented in 
[13]. As to the attractivity of the FP, we have found that, 
when it exists, there are always finitely many attractive 
directions. Therefore gravity seems to remain asymptot- 
ically safe also in the presence of generic matter fields. 
We expect that this result will still hold if we add other 
interactions between matter fields that are asymptoti- 
cally free. From this point of view, the scalar field posed 
a greater challenge, since the pure scalar theory is not 
asymptotically free. It is remarkable that the coupling 
to gravity fixes this problem and at the same time also 
offers a solution to the triviality problem. 

All these results add on to the other proofs that have 
been collected in the literature about the physical relia- 
bility of this approach. 
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